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(54) Method for determining 
distance and direction from an open 
well to a cased well 

(57) A method for determining 
distance and direction from an open 
well to a cased well comprises: (1) 
utilizing resistivity logs obtained in the 
open well, and/or in the cased well 
before casing in run, to determine 
characteristics of earth strata Z to Zn 
penetrated by the wells, (2) 
determining the actual trajectory of 
the open well from survey data, (3) 
determining theoretically expected 
resistivity curves for long and short 
spaced electrode-type tools A, M, N at 
various depths in the open well for 
different assumed positions of the 



cased well relative to the open well, 
(4) determining measured resistivity 
curves for tools A, M, N by logging the 
open well, (5) varying the assumed 
relative positions of the open and 
cased well to adjust the spacing of the 
curves until the expected and 
measured curves match in a best fit 
mathematical sense; the amount of 
variation in assumed relative position 
required to match the curves 
determining the distance of each 
station in the open well to the cased 
well, and (6) determining the direction 
of the cased well from the open well 
by constructing a plan view of the 
open well trajectory from the survey 
data and utilizing the distances 
previously determined. 
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SPECIFICATION 

Method for determining distance and direction from an open well to a cased well 

This invention relates to a method for utilizing well log data and other measurements to 
determine the distance and direction between a cased blown out well and an adjacent open (e.g. 
5 uncased well) relief well. The intersection of a blown out well by a relief well is frequently required in 5 
order to kill the wild well by pumping cement or other fluids into its well bore. It is an extremely difficult 
and complex engineering problem to effect an intersection in the subsurface at depths frequently 
greater than 1 5,000 feet. This problem is much simpler if a method is available for determining the 
distance of the relief well to the casing during the course of its approach to the target well. The 
1 o invention is also concerned with the use of similar data and measurements to avoid having two wells 1 0 
intersect, which can be a problem where the wells are drilled from a common starting point such as an 
off-shore drilling platform. 

A number of methods have been proposed for utilizing well logging data to estimate distance and 
direction from an open well to a cased well. Direction and distance to casing can be estimated by 
1 5 lowering a magnetometer on a wireline into the relief well and measuring the magnetic field and/or its 1 5 
gradient produced by the natural magnetization of the casing. This method is described in U.S. Patent 
No. 4,072,200. 

A similar method is described by J. D. Robinson and J. P. Vogiatziz in U.S. Patent No. 3,725,777 
and in a publication entitled "Magnetostatic Methods for Estimating Distance and Direction from a 
20 Relief Well to a Cased Borehole", Journal of Petroleum Technology, v. 24, no. 6, pp. 741 — 750, June 20 
1972. 

A limitation of these methods is that the magnetic field and/or its gradient can be too weak for 
accurate detection unless the magnetometer is near the end of the casing where the field is strongest. 
This method is a passive method since it depends on the casing having a natural magnetization. 
25 Professor Arthur F. Kuckes of the School of Applied and Engineering Physics at Cornell University has 25 
recently developed two active methods for estimating distance and direction from a relief well to a 
cased target well. 

The first method developed by Professor Kuckes and co-workers involves the creation of a low 
frequency alternating current in the casing of the blown out well by placing electrodes on the surface in 

30 the ground near the wild well. The conductivity of the metal casing is roughly 1 0 6 times greater than 30 
that of the earth so that the current finds it favourable to flow in the casing. Thus, having established an 
alternating current in the casing the magnetic field produced by this current is measured by a 
magnetometer placed in the relief well on a wireline cable. Suitable interpretation of the direction and 
magnitude of the measured magnetic field can be used to estimate both distance and direction to the 

35 cased target well. A practical limitation of this method occurs because a blow out frequently results in 35 
the separation of casing joints. A result of this is that current leakage from the casing can prevent the 
induce^ currents from reaching desired depths with sufficient strength to produce measurable and 
interpretable magnetic fields. 

Professor Kuckes and co-workers have also developed a second active method. This method 

40 places both the source current electrodes and the magnetometer in the relief well on a wireline logging 40 
cable. In this case the magnetometer measures a reflected magnetic field due to the presence of the 
casing. This method has serious interpretation problems if the relief well trajectory is not a straight line 
path. In practice, it is more common for a well to have a tortuous trajectory in space over distances of 
hundred to thousands of feet. In this case the electrodes producing the current and the magnetometer 

45 will not be on the same straight line. This results in a primary and secondary magnetic field at the 45 
magnetomer. The primary field would be zero if the relief well trajectory were a straight line path. 

In meandering wellbores the total magnetic field at the magnetometer consists of both the 
primary field and the secondary field. The secondary field results from the presence of the casing and 
therefore it is this field which contains information about the djstance and direction of the relief well to 

50 the blown out well. In relief wells with tortuous paths it may not be possible to separate the primary 50 
and secondary fields with sufficient accuracy to obtain accurate distance and direction estimates using 
this method. 

A different method which is more closely related to the method of this invention is disclosed in 
U.S. Patent No. 3,748,574 and described in a publication by F. R. Mitchell et al entitled "Using 

55 Resistivity Measurements to Determine Distance Between Wells" printed in the Journal of Petroleum 55 
Technology, v. 24, no. 6, pp. 723—741 , June, 1 972. This patent and publication represent the state of 
the prior art most relevant to the present invention. The method of the Mitchell et al publication utilizes 
logging data obtained in the relief well with an Ultra-Long Spaced Electric Log (ULSEL), electrical 
survey (ES) and induction and/or other shallow investigating resistivity logs. The latter logs are 

60 common and familiar to operators engaged in the drilling of oil and/or gas wells. The ULSEL tools were 60 
developed and are offered as a service of Schlumberger Well Surveying Corporation. These tools are 
electrode-type electrical logging tools which have electrode spacings of hundreds to thousands of feet. 
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A method for using the ULSEL data to detect salt domes at distances of hundreds to thousands of 
feet from a well has been described in U.S. Patent No. 3,256,480. This patent contains a description of 
the ULSEL tools. The application of the ULSEL tools to salt dome detection was the original impetus for 
the development of the tool. Salt domes, as is well known, represent structural traps for hydrocarbons 
5 and therefore their presence, location and size are of obvious interest to geologists. The ULSEL can 5 
detect the presence of salt domes because these structures have anomalously high resistivities as 
compared to other geological structures present in the subsurface. 

The aforementioned publication by Mitchell et al and U.S. Patent No. 3,748,574 describe a 
method and field examples wherein ULSEL, ES and induction log data are utilized to estimate distance 
10 of a relief well to a cased target well. The ULSEL can be used to detect the presence of casing since in 1 0 
proximity to casing the resistivities measured by the ULSEL are reduced significantly relative to their 
values in the absence of casing. This reduction in measured resistivities when the casing is within the 
range of investigation of the ULSEL is a result of the fact that the resistivity of the casing is roughly 
1 0" 6 times that of the earth. 

1 5 Mitchell et al discuss a method for estimating distance to casing by utilizing the ratio of measurec 1 5 

to expected (e.g. theoretically computed in the absence of casing) resistivity. This ratio decreases as 
the casing is approached by the relief well. A method is proposed for using the aforementioned 
resistivity ratio to estimate distance to the cased target well. The method proposed by the 
aforementioned Mitchell et al publication and U.S. Patent No. 3,748,574 can be described as follows. 

20 The relief well is logged with an induction tool in an interval of interest which in practice could range 20 
from tens to thousands of interest which in practice could range from tens to thousands of feet. Since 
the depth of investigation of induction logging tools is relatively shallow, the measured resistivities are 
not affected by proximity to casing provided that the relief well is greater than roughly twelve feet from 
the target well. The induction log resistivities are used to make a layered resistivity model of the earth 

25 which best approximates the actual resistivity versus depth profile determined from the induction log. 25 
The number of layers included in the model is variable and in practice could be several hundred. 

The effects of resistivity anisotropy are included in the layered model by having, in each layer, a 
resistivity parallel to the layering planes and a resistivity perpendicular to the layering planes which are 
not equal in an anisotropic layered earth. Since only a finite number of layers can be included in the 

30 model, it is assumed that the layered medium is bounded above and below by semi-infinite half spaces 30 
with appropriately chosen resistivities. The theoretically expected response of the ULSEL is calculated 
for this layered medium. This then gives the theoretical response of the ULSEL in the absence of the 
cased target well. The method of Mitchell et al assumes that the relief well trajectory is a straight line 
and therefore that the direction of the relief well has a constant bearing (e.g. angle measured relative tc 

35 north) and deviation (e.g. angle measured relative to the vertical). A limitation of this method is that in 35 
practice a well borehole frequently follows a meandering trajectory in space. Failure to take into 
account this trajectory leads to errors in the theoretically computed ULSEL responses and, therefore, to 
errors in distances estimated from the ratio of measured to expected resistivities. ... 

Another limitation of this method is that directional information cannot be obtained using the 

40 assumption of a straight line path for the relief well. The second step in the Mitchell et al method . 40 
utilizes the ratio of measured to expected resistivities to estimate distance of the relief well to the 
cased target well. This step of the method makes an unrealistic and inconsistent approximation. The 
approximation is that the cased target well is not in a layered medium but is situated instead in an 
infinite homogeneous medium. This approximation is not valid, is inconsistent with the first step of the 

45 method which constructs the layered model, and therefore it introduces unknown errors in the 45 
distances estimated using this method. 

It is an object of the present invention to avoid or minimize one or more of the above 
disadvantages, using ULSEL, ES, induction and/or other shallow investigating resistivity logs and 
directional survey data and measurements to account for the trajectory of the relief well and the fact 

50 that the casing is in a layered earth in determining distances frorn {he relief well to the cased target 50 
well. Moreover, for a tortuous relief well trajectory, the method disclosed herein provides means for 
determining direction from the relief well to the target well in addition to distance. 

A method has now been discovered for utilizing resistivity and directional survey log data in order 
to determine more accurately than heretofore possible distances from a relief well to a cased target 

55 well. The method takes into account the actual trajectory of the relief well in space and therefore 55 
provides directions as well as distances to the target well. 

The present invention provides a method for determining distance and direction from an open 
well to a cased well comprising determining the resistivities and thicknesses of the earth strata 
penetrated by said wells by induction logs run in said open well or previously run in said cased well 

60 determining the actual trajectory in the earth of the open well by directionally surveying said well, 60 
determining theoretically expected resistivity curves for long and short spaced electrode-type resistivity 
tools at various depths in said open well for different assumed positions of the cased well relative to 
the open well, from said earth strata resistivities and thicknesses determined by induction logs, 
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translation or saia theoretically expected resistivity curves until the expected resistivity curve matches 
the measured resistivity curve in a least squares fit mathematical sense, and determining the distances 
of selected points in the open well to the cased well from the amount of translation of said theoretically 
expected resistivity curve required to produce said match of said theoretically expected and said 
5 measured resistivity curves. 5 
The increased accuracy of this method results from essentially two improvements over prior art 
methods. Firstly in making this invention, it has been discovered that greater accuracy in distance 
determinations is obtained by incorporating into the method means for taking into account the actual 
trajectory of the relief well in space. The method determines, at each station depth (e.g. chosen at one 
1 o foot intervals) in the relief well, the spatial co-ordinates of the stations over a specified depth interval. 1 0 
Thus the actual positions of the ULSEL electrodes "A", "B", "M" and "N" in the relief well are known 
accurately. This knowledge is essential in order to determine accurate values of expected resistivities 
which are then compared with measured resistivities to obtain distances and directions to the target 
well. The assumption of a stright line path, in prior art methods, precluded using resistivity 
1 5 measurement to obtain any directional information. 1 5 

The method of this invention also provides greater accuracy in determining distances than prior 
art resistivity methods because the method avoids the inconsistency and inaccuracy present in the 
prior art. This inconsistency and inaccuracy in the prior art is a result of the procedure for determining 
distances to the target well. The prior art procedure assumes that the cased target well is in a 
20 homogeneous medium rather than in the layered earth. This assumption is invalid and leads to 20 
unknown errors in distances determined from prior art resistivity methods. The method of this 
invention employs a procedure for determining distances which automatically takes into account the 
fact that the cased target well resides in a layered medium where the resistivities of the layers are 
determined from the induction log and/or other shallow investigating resistivity logs. The distances to 
25 the target well are determined by varying the assumed position of the target well relative to the relief 25 
well in order to shift the relative positions of the measured and the expected resistivities. The correct 
relative position is obtained when the expected (e.g. theoretically calculated for the cased target well in 
the layered earth) ULSEL and ES resistivities match the measured ULSEL and ES resistivities. The 
distances to casing of all stations in the relief well are then accurately determined. 
30 Further features and advantages of the invention will become more readily apparent from the 30 

following detailed description given by way of example and illustrated with reference to the 
accompanying drawings. 

Figure 1 is a schematic cross-section showing a cased and open well in a layered earth; 
Figure 2 is a schematic cross-section showing an ULSEL tool with electrodes "A", "M" and "N" 
35 in an open well; 35 
Figure 3 illustrates the relation between the expected and measured ULSEL (1 50/350) 
resistivities for different assumed positions of a cased well relative to an open well; and 

Figure 4 is a schamatic plan view of an open well trajectory illustrating means for determining the 
direction of a cased well relative to the open well. 
40 In Figure 1 is shown a schematic cross-section of the earth showing a portion of a relief well path 40 

as the well approaches a cased target well for a planned intersection. The first step in the method 
requires determination of the trajectory of the relief well over the depth interval of interest. The drill 
string is removed from the relief well and a borehole directional survey is taken at various points P,, 
P 2 , . . .,P n in the relief well. At each point the directional survey data consists of a measured depth, a 
45 bearing angle and a deviation angle. These two angles specify the direction and orientation of the relief 45 
well trajectory at the measured depths. 

It would be desirable to have directional survey data at one foot intervals over the depth interval 
of interest (e.g. from P 1 to P n ). In practice, the directional survey data is generally not measured at such 
closely spaced intervals. The direction of the relief well at one foot intervals along its path from P 1 to 
50 P n , however, can be obtained by interpolating between the measured directional survey data points. 50 
Interpolation techniques have been developed and are discussed^ for example, in a publication entitled 
"Directional Survey Calculation" by J. T. Craig, Jr. and B. V. Randall in Petroleum Engineer, March 
1976. 

An induction log is run in the relief well either before or after the directional survey. The induction 
55 log resistivity data is used to construct a layered model for the earth. The apparent thicknesses along 55 
the path and the resistivities for the respective layers are chosen to best approximate the resistivity 
profile obtained with the induction tool. Previous induction log rungs in the relief well are used to 
obtain thicknesses and resistivities for depths above P,. The layering structure depths below P n can be 
obtained from induction logs run in the target well. Generally the layering structures above P, and 
60 below P n do not have a strong influence on the directions and distances to casing determined by this 60 
method. In practice it will frequently suffice to approximate these upper and lower zones by infinitely 
thick layers with suitably chosen resistivities. In each layer the effects of resistivity anisotropy are taken 
into account by allowing for different formation resistivities in the directions parallel and perpendicular 
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from P, to P n , the cartesian co-ordinates (x'„ y'„ z',) of each station are determined relative to an origin 
0' chosen at one of the stations, say at station P,. The index i ranges from 1 to n s where n s is the ^ 
number of stations chosen at one foot intervals between P, and P n . 

In Figure 2 is shown a schematic cross-section of the earth showing an ULSEL or ES topi situated 
in the relief well. The tools consists of a current electrode J "A" and two voltage measuring electrodes 
"M" and "N". The current return electrode "B" is not shown in Figure 2 since, in practice, it is 
frequently located on the earth's surface in the mud pit. In this case one can assume, for practical 
purposes, that the spacing AB between the A and B electrodes is infinite since the B electrode does not 
influence the expected (e.g. theoretically determined) resistivities. 

For the short spaced ES tools such as the 1 6-inch and 64-inch Schlumberger Normal devices the 
B electrode is in the borehole and has an effect on the expected resistivities. It is well-known to those 
skilled in the art of electrical logging that it is straightforward to include the effect of the B electrode in 
determining expected resistivities. Since the inclusion of a B electrode in the borehole would make 
certain formulas given below, unnecessarily complex and is not relevant to the present invention, the 
discussion given below in conjunction with Figure 2 will, for simplicity, assume that the B electrode is 
at infinity. In Table 1 are listed some of the common electrode spacings for ULSEL and ES tools. 

Table 1 

ULSEL and ES electrode spacings 



20 



25 



30 



Tool 


Electrode spacings (feet) 


ES Tool 


AM 


AN 


AB 


1 6-inch Normal 


1.33 


20 


89 


64-inch Normal 


5.33 


71 


89 


20-foot Normal 


20 


71 


infinity 


ULSEL Tool 








75/350 


75 


350 


infinity 


75/600 


75 


600 


infinity 


150/350 


150 


350 


infinity 


1 50/600 


150 


600 


infinity 


,150/1000 


150 


1000 


infinity 


600/2400 


600 


2400 


infinity 


1000/2400 


1000 


2400 


infinity 
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The method of this invention for determining distances and directions from the relief well to the 
cased target well involves moving the assumed location of the casing relative to the relief well 
trajectory which has been determined, as discussed above, using the directional survey data. For each 
specified location of the target relative to the relief well, the expected resistivies of the various ULSEL 
and ES devices are determines. The correct location of the target is obtained when the expected 
resistivites match the measured resistivities obtained by logging the relief well with the same ULSEL 
and ES devices for which the expected resistivities have been determined. 

The following theoretical consideration will make clearer the procedural steps of the method. 
First the expected resistivities are determined for a given (e.g. specified) position of the target well 
relative to the relief well. Referring to Figure 2 a schematic is shown of an ULSEL or ES tool having 
electrodes A, M and N situated in the borehole of the relief well. Also shown in Figure 2 is the origin 0 
of a co-ordinate system fixed in the target well. Recall that the cp-ordinates(x' if yVz'j) of the relief well 
are known from the directional survey data at one foot intervals from the points P, to P n in Figure 2. 
These primed station co-ordinates are relative to an origin 0' fixed at one of the stations on the relief 
well path. — 

Let the origin 0' have cartesian co-ordinates x 0 , 



relative tc the origin O fixed in the cased 



target well. The co-ordinates x 0 , y 0 . z 0 . describe the location and position of the target well relative to 
a •• -'of each station in the relief well relative to the origin 



the relief well. The cartesian co-ordinates (x„ y 
0 fixed in the casing are: 



x,=x',+x 0 



v,=y',+y<,. 



(1a) 
(1b) 



35 



40 



45 



50 



and 
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where the index i=1 , 2, . . . n s is equal to the number of stations situated at one foot intervals along the 
relief well path. The equations or formulae used throughout this description are numbered 
parenthetically for purposes of reference. 

The distances of each station to the target well are given by 

5 dHix'i+xJ^ly'i+y,,) 2 ) 0 ' 5 (2) 5 

The method of the present invention determines accurate values for (x 0 , y 0 ) and therefore for the 
distances from the relief well stations to the cased target well. The method involves "moving" the 
target well's location theoretically relative to the relief well trajectory by varying the choice for the co- 
ordinates (x 0 , y 0 ). The correct values of (x 0 , y 0 ) are found when the expected resistivities determined for 
1 0 the various ULSEL and ES tools match the measured resistivities obtained from logging these tools in 1 0 
the relief well. 

The correct values for (x Q , y 0 ) having thus been determined, the distances d, of each station in the 
relief well are obtained from formula (2). The greater accuracy provided by the present method over 
prior art methods is due, in part, to the fact that this method determines expected resistivities for the 

1 5 ULSEL and ES tools, for different assumed positions of the cased target well located as shown in Figure 1 5 
2 in a layered medium. The effect of both the casing and the layering on the expected resistivities is 
therefore properly accounted for by the present method. The requirement that the expected ULSEL and 
ES resistivities match (in a best fit sense described below) the measured values obtained by logging the 
relief well then results in more accurate determinations of distance to the target well than are obtained 

20 using the prior art method. 20 
In order to distinguish the present method from the prior art method discussed in U.S. Patent No. 
3,748,574 it is necessary to recall the prior art method. The method of patent 3,748,574 determines 
expected resistivities for the ULSEL and ES tools in a layered medium without taking into account the 
effect of a cased target well on the expected resistivities. That is, the layered medium for which the 

25 expected resistivities are determined does not contain a cased target well. The proximity to the cased 25 
target well is then inferred when the measured resistivities in the relief well are less than the expected 
resistivities. The determination of distance from the target well to the relief well is made by assuming 
that the target well is located in a homogeneous (i.e. nonlayered) medium. Thus, that method is self- 
contradictory and the distances determined therefrom are subject to unknown uncertainties. 

30 Implementation of the method of this invention requires the ability to determine the expected 30 

resistivities for ULSEL and ES tools in a layered medium containing a cased target well. In making this 
invention, the following formulas were derived for this purpose. The expected resistivities (e.g. also 
known as apparent resistivities) are determined by solving the Poisson's equation obeyed by the 
electrostatic potential <p (p, 0, z) where (p, 0, z) are cylindrical co-ordinates relative to an origin 0 (e.g. 

35 see Figure 2) fixed in the cased target well. 35 
Tbe formulation of the theory of expected (e.g. apparent) resistivities is well-known and, for 
example, is discussed by V. N. Dakhnov in Geophysical Well Logging, Moscow Petroleum Institute 
( 1 959), translated by G. V. Keller in Colorado School of Mines Quarterly, Vol. 57, No. 2 ( 1 962). Using 
the above referenced work and a paper by D. H. Shun entitled "Green's Functions for Composite 

40 Media", in the International Journal of Engineering Science, Vol. 16, pp.475 — 482 (1978), the 40 
formulae given below can be developed by those skilled in the art of theoretical physics. 

The electrostatic potential <p(p, 0, z) at any point in the layered medium containing the relief well 
and the cased target well can be obtained from the formula, 

oo ^ 

$(p e f z)-/d\A2 ^cos n(e-e) F n ( \ ) Gj < n > (a), o) 



45 where £ n =1 for n=o and c n =2 for n> 1 . 

In formula (3) the functions F n (A) and G, ,n, (A) are defined by the formulae 



n ,2 ■ v/* 



Jf,(Aa) + Y'Ua) 



and 



45 
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where. 



Aj (n) (A)--i— ^ } (J n (p 0 a ) Y n (xa) - Y n (Ap 0 )j n (Ad) } 



(6) 



and. 



J (z Q )= 2Aa m (b m c m -a m d m ) . 



(7) 



5 The functions J n and Y n in formula (4) are ordinary Bessel functions of order n of the first and 5 

second kinds, respectively. In the above formulae (/> 0 - f 0 O ' z o> are the cylindrical co-ordinates of the 
current electrode A in the relief well, I is the strength of the current emanating from A and "a" is the 
outer radius of the target well casing. The functions Gj ,n, (A) are defined for Z;<z<z i+1 (e.g. z in the ith 
layer) with i=1 , 2 . . .N where N is the number of layers as shown in Figure 2. In formula (5) the 
1 0 subscript m denotes the layer in which the A electrode is located (e.g. z m <z 0 <z m+1 ). The coefficients a r 1 0 
bj, Cj and dj are determined from boundary conditions which require continuity of the potential and the 
component of the current density perpendicular to the layers at each interface. The 2 (N— 1 ) 
coefficients a t and bj are obtained, for example, from the following recursion formulae, 



'b.. 



15 



and 



(i) (aj+aj_|)Az 



(i) (dj-aj 



i-l 



I — 



(8) 



(9) 



15 



for i=2, 3 N, where a, is arbitrary and b^O. 

The coefficients M, 01 and M 2 (i) are defined as follows: 



M ''11/2 0*°'-'^''') 
1 a, (SyT 1 ' 



20 



and 



m 2 -1/2 (i- '-' V . . . ) 



(IO) 



(ii ) 



20 



where 



(0 



(i) rt o-s 



(12) 



with crJ" and a 111 the conductivity in the horizontal direction in the ith layer and the conductivity in the 
25 vertical direction, respectively. In the above formulae the z, are the co-ordinates of the interfaces of the 2 5 
layers. The coefficients c, and d, obey backward recursion formulae and are given by: 



c,_, - M 3 (l) « (a r a i-l )Az i c .+ 



14 



(i) - (a " +a J-l ,A2 i d 



i » 



for i=N, N— 1 , ...2, with c„=0 and d„ arbitrary and where 



03! 



(14) 



30 



M 



(i) 



and 



l/2 (I 



l/2 (I- 



(i) 



v \ 
TT=T) ) 



(I5) 



30 



(i) 

I 

Trrrr 



(i6) 
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The above formulae can be used by those skilled in the art of computer programming and 
numerical analysis to determine expected (e.g. apparent) resistivities for ULSEL and ES tools having 
various electrode spacings. The expected resistivity R a at the location of electrode A is determined from 
the formula, 

5 ^^Wl^tlVD-^NIKI/UAMl-l/UANir 1 (17) 5 

where <p(M) and <p(N) are the respective potentials determined from the above formulae at electrodes M 
and N. The quantity L(AM) is the straight-line distance (in general not equal to the electrode spacing) 
separating the A and M electrodes in the relief well. Likewise, the quantity L(AN) is the straight-line 
distance (in general not equal to the electrode spacing) separating the A and N electrodes in the relief 

1 0 well. The distances L(AM) and L(AN) are known accurately since the co-ordinates of each station in the 1 0 
relief well are determined from the directional survey data (or by its interpolation). 

The expected resistivities are computed for different locations of the A electrode in the relief well 
for ULSEL and ES tools having different electrode spacings and for different assumed positions (e.g. 
different values of (x 0 , y 0 » of the target well relative to the relief well. The computations of the expected 

1 5 resistivities taken into account: ( 1 ) the actual path of the relief well in space as determined from the 1 5 
directional survey data and therefore the actual locations in three-dimensional space of the A, M and N 
electrodes, (2) the effect of the cased target well situated in the layered medium containing the relief 
well for different assumed relative positions of the two wells, (3) the effect of resistivity anisotropy for 
the different layers in which the two wells are located, (4) the effects of different electrode spacings 

20 and configurations and (5) the effects of different target casing radii. 20 
The expected resistivites are determined for different assumed positions of the target well relative 
to the relief well and plotted versus the A electrode measured depths in the well. The three cases 
illustrated in Figure 3 show the effects of varying the assumed position of the relief well relative to the 
target well. For case 1 the position assumed for the target well is closer than the actual position of the 

25 relief well to the target well. The expected resistivities are therefore less than the measured 25 
resistivities. For case 2 the assumed position of the target well relative to the relief well is farther than 
the actual relative positions. Therefore the expected resistivities are greater than the measured 
resistivities. For case 3 the assumed position coincides with the actual relative position and the 
measured and expected resistivities match. A quantitative measure of the goodness of the match can 

30 be obtained by minimizing a function S k defined by 30 
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mk 



where k denotes a particular set of ULSEL or ES spacings R aK li) and R mk ,M are the expected and 
measured resistivites, respectively, for the A electrode at station i in the relief well. The values of (x Q , y Q ) 
which minimize S k determine the distances from each station in the relief well to the target well. In 

35 practice ULSEL and ES tools having different electrode spacings are utilized. The longer spaced tools 35 
are used to provide a long range determination of distance to casing whereas the shorter spaced tools 
are used when the target is within their range of investigation. 

In Figure 4 is shown a schematic plan view of the relief well trajectory as determined from the 
directional survey data. The direction to the target well is determined for a tortuous relief well path as 

40 shown in Figure 4. Direction to the target well is determined by choosing at least three points along the 49 
relief well path. Let these points be centres of circles whose radii are the distances to the target well 
which have been determined as described above. The common intersection of these circles defines the 
direction of the target well from the relief well. It is necessary t{iat the relief well p&th be curved in 
order that this method determine the direction uniquely. In practice; it is advantageous to.intentionally 

45 drill a curved path as the relief well approaches the target in order to make the determination of 45 
direction more accurate. 

The following example illustrates the application of this method to a specific case. A cased well is 
blown out and flowing from a zone at 1 3,500'. The casing in the wild well is set to a depth of 1 6,045'. 
A directional well is drilled with the objective of intersecting the wild well at 1 3,500' in order to kill the 

50 flow by pumping cement and other fluids into the flowing zone. A string of protective casing is set in 50 
the relief well to a depth of 1 1 ,250'. 

The approach to the intersection is done in steps. At the first step, the relief well total depth (TD) 
is 1 2,650'. The protective casing in the relief well can adversely perturb the ULSEL measured 
resistivities. The interval where this occurs is roughly an AM+AN spacing below the casing depth. 

55 Therefore, the interval where the measured ULSEL resistivities are valid extends from an AM+AN 55 
spacing below the protective casing to the TD of the relief well. 
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being cased, and induction log data obtained in the wild well below 1 2,650' prior to its being cased are 
utilized to construct a layered resistivity model of the earth strata penetrated by the wells. The layered 
model assigns resistivities and thicknesses to the earth strata. 

The trajectory of the relief well in the interval from 1 1 ,250' to 1 2,650' is determined from a 
5 borehole directional survey. The survey data are taken at 50' intervals in the relief well and are 5 
interpolated using a computer programme to obtain the spatial co-ordinates (relative to an origin in the 
relief well) of the relief well at one foot intervals from 1 1 ,250' to 1 2,650'. 

The layered resistivity model and relief well co-ordinates are used as input into formulae (3) — 
( 1 7) and formulae (1 ) to compute the expected ULSEL and ES resistivities. The computations are done 

10 forelectrodespacings of 20/7 1,75/150, 75/350,75/600, 150/350 and 1 50/600. For each set of 10 
spacings, the expected resistivities are calculated for different positions of the A-electrode in the relief 
well and for different assumed relative positions (e.g. different choices for (x D , y Q ) of the relief and wild 
well). The computations are done by a Control Data Corp. (CDC) CYBER 1 76 computer and the 
expected resistivities are plotted versus depth using software available on the CDC computer. 

1 5 First, it is determined, for each of the ULSEL and ES tools, whether the wild well is within the 1 5 

range of investigation of the tool. This is done by first assuming that (x Q , y 0 ) are at infinity. The expected 
resistivities versus depth are then computed for each of the tools and compared with the measured 
resistivities obtained by logging the tools in the relief well. In the present example, the measured and 
expected resistivities show no consistent differences for the 20/71 and 75/1 50 tools. Thus, it is 

20 concluded the wild well is outside the range of investigation for those two tools. For the 75/350, 20 
75/600, 1 50/350 and 1 50/600 spacings, the measured resistivities are consistently lower than the 
corresponding expected resistivities which implies that the wild well is within the range of investigation 
of each of these tools. 

The computations are repeated for the 75/350, 75/600, 1 50/350 and 1 50/600 tools using 

25 different choices for <x 0 , y 0 ). For example, in case 1 of Figure 3, there are plotted the measured and 25 
expected resistivities for the 1 50/350 tool for assumed values of x 0 =60' and y 0 =30'. The expected 
resistivities are consistently less than the measured values, indicating that the actual position of the 
wild well is farther from the relief well than the assumed position, in case 2 of Figure 3, results are 
shown assuming that x o =200' and y 0 =80'. Note that the expected resistivities are consistently greater 

30 than the measured values, indicating that the wild well is closer to the relief well than assumed. 30 
In case 3 of Figure 3, the measured and expected values are plottea for x o =90' and y o =60'. There 
are no consistent differences in measured and expected resistivities and therefore the assumed relative 
position is correct. Using the values x o =90' and y o =60', the distance from each station in the relief well 
to the wild well can be obtained from formula (2) and the interpolated directional survey data. For 

35 example, the first station irj the relief well which is taken as the origin 0' (e.g. x' t =y',=0) is a distance of 35 
1 08' from the wild well. Confirmation of the correct assumed relative position is obtained by 
computing the expected resistivities for the 75/350, 75/600 and 1 50/600 tools using x 0 =90' and 
y o =60'.The absence of consistent differences between the expected and corresponding measured 
resistivities supports the results obtained with the 1 50/350 tool. 

40 The directional survey data are utilized to construct a plan view of the relief well trajectory. A? 40 

illustrated in Figure 4, the direction to the relief well is determined by choosing at least three points on 
the path as centres of circles whose radii are the distance of the respective centres from the wild well. 
The common intersection of these circles gives the direction of the wild well from the relief well. 

The directional information is utilized to direct the relief well toward the wild well as the drilling of 

45 the relief well is commenced. After drilling several hundred feet further, the above steps are repeated 45 
and necessary changes made in the direction of the relief well to direct it to intersection with the wild 
well. 

The measurements and calculations made as described above can be used in a like manner for 
the purpose of keeping two wells from intersecting, as in t£e case of wells drilled from closely spaced 
50 positions such as an off-shore drilling platform. .< 50 

While this invention has been described fully and completely - With special emphasis upon a single 
preferred embodiment, it should be understood that within the scope of the appended claims, the 
invention may be practiced otherwise than as specifically described herein. 

Claims . . 

55 1 . A method for determining distance and direction from an open well to a cased well comprising 55 

determining the resistivities and thicknesses of the earth strata penetrated by said wells by induction 
logs run in said open well or previously run in said cased well determining the actual trajectory in the^ 
earth of the open well by directionally surveying said well, determining theoretically expected resistivity 
curves for long and short spaced electrode-type resistivity tools at various depths in said open well for 

60 different assumed positions of the cased well relative to the open well, from said earth strata 60 
resistivities and thicknesses determined by induction logs, determining measured resistivity curves for 
long and short spaced electrode-type resistivity tools by logging said open well, varying the assumed 

-„l^+;w^ rt «*»;+;^r»e nf tho nnon anH racori \a/oII hv tranclatinn of <saiH thpnrPtinallv ftXHRCtftd resistivity 
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mathematical sense, and determining the distances of selected points in the open well to the cased 
well from the amount of translation of said theoretically expected resistivity curve required to produce 
said match of said theoretically expected and said measured resistivity curves. 

2 A method according to claim 1 including determining the direction of said cased well from said 

5 ooen well by constructing a plan view of the open well trajectory from said directional survey data, 5 
selecting at least three points on the trajectory as centres of circles with the radius of each circle equal 
to the previously determined distance of its centre from the cased well, and determining the location of 
the cased well relative to the open well by the common intersection of said circles. 

3 A method according to claim 1 or claim 2 wherein said spaced electrode type resistivity tools 

1 o include a current electrode A and two voltage measuring electrodes M and N spaced longitudinally in i o 
said open well, the expected resistivity at the location of electrode A in the open well being determined 
by the formula: 

R a =4 Tr/MpdvD-^N)) ( 1/UAM)-1/UAN))-' 

where I is the current emanating from the A electrode, L(AM) and L(AN) are the straight line distances 
1 5 from the A electrode to the M electrode and the N electrode, respectively, as determined from the 1 5 

directional survey data, <p(M) and <p(U) are the respective electrostatic potentials at the locations of 
electrodes M and N, the electrostatic potential cpip, 0, z) at any point {p, 0, z) in the layered earth 
containing both the open and cased well is determined by the formula: 

oo 

+ (p e, z)-/dJUl ^cos n(e-e) F n ( A ) Gj < n ) (a), 

o 

20 where £ n =1 for n=o and £ n =2 for n > 1 , functions F n '' l > and G^'IA) are defined by the formulae: 20 

j n Up)Y n U a)-J n (A a)Y n (a p) 



F n (A) - 

j£(Aa) + Y n 2 Ua) 



n TTTTTTZT 



and 




where, 

(n)m I 



25 A i U,= "2TTlm) (J n(PoA> Y n<Aa)-Y n (A Po )j n (Aa)) 



25 



and 

J (z Q )= 2\a m (b m c m -a m d m ) 



where J and Y n are ordinary Bessel functions of order n of the first and second kinds, respectively, 
where (o 6 a . zj are the cylindrical co-ordinates of the current .ejectrode A in the relief well. I is the 
30 strength of the current emanating from A and "a" is the outer radiusof the target we I casing, where 30 
the functions G' n, (A) are defined for z in the ith layer so that z,<z<z l+1 with i=1 , 2 . . . N, where N is the 
number of layers being measured, the subscript m denotes the layer in which the A electrode is 
located, z^^Zo^z^,, coefficients a„ b„ c, and d, are determined from the following recursion formulae, 

o, -M l (l ) 4 - (a i- a '-' , %-,^ il e- ( ^ a .-l , "i bj . | 



35 and 



b f -M 2 ,i, c ,a i ta -' , "V 1 ^' i, «' arQi - |UZi '> i -, 
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Ml ,i '-,/2,u2i=L2s^l , ) 

a,- <f v 0 ) 

and 

M 2 -l/2 (I- V ■ !■ ) 

' a,a v (i) 

where 

with cr h ( " the conductivity in the horizontal direction in the ith layer and a v m the conductivity in the 
vertical direction* in the above formulae the z, are the co-ordinates of the interfaces of the layers, 
coefficients c f and d, are determined by the formulae: 

1 ° and , „ w , 10 

d|-, .M 4 ' , >«^* 0 W>» I lq + M3<"c" (arC, - |lM 'di. 

for i=N, N-1 , . . .2, with c N =0 and d N arbitrary and where 

(i) ,. <>i «v (il 



and 



15 



a. <J <') 



4. A method according to any one of claims 1 to 3 wherein the expected resistivity curves are 
determined by calculations based on the actual trajectory of the open well in the earth as determined 
by interpolation of directional survey data taken in the well. 

5. A method according to any one of claims 1 to 3 wherein the expected resistivity curves are 

20 determined by calculations involving the present of a "B" electrode in the borehole of the open well. 20 

6. A method according to any one of claims 1 to 3 wherein the expected resistivity curves are 
determined for different assumed positions of the cased well relative to the open well by. varying 
cartesian co-ordinates (x 0 , y Q , z Q ) which describe the relative positions of the cased and open wells. 

7. A method according to claim 6 wherein cylindrical co-ordinates (p x , 6,, z,) of stations in the 

25 open well relative to an origin 0 in the cased well are defined by the formulae: 25 

P,-((x; + x Q ) 2 + Gy', - y 0 ) 2 ) 0,5 



1 

and 



e ; - tan" 1 t y 'i + y o } 
( x'j ■»- x Q ) 



Zj - z'j + z Q . 



30 for i=1 , 2 n s , where n s is equal to the number of stations at one foot (30 cm) intervals in the open 30 

well and where the co-ordinates, x ; ,, v',, z' ; , are determined by interpolation of the directional survey 
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8. A method according to claim 6 or claim 7 wherein the expected resistivity curves are 
determined for short and long spaced electrode-type resistivity tools having selected spacings and 
where the assumed values (x 0 , y 0 , z 0 ) are varied in order to minimize the function S k given by the 
formula: 



S k - I — 211 

mk 



i-l (Rl»J) 2 



where the index k denotes a particular set of electrode spacings, rj" and R mk « are the expected and 
measured resistivities, respectively, for the A electrode at the ith station in the open well. 

9. A method according to any one of claims 6 to 8 wherein the particular values of (x 0 , y e , z 0 ) 
which minimize S k determine the distances to the cased well of each station i in the open well 
1 0 according to the formula: 

d l =((x',+x 0 ) 2 +(y',+y.) 2 )" 



1 0 A method according to any one of claims 3 to 9 including determining the direction of said 
cased well from said open well by constructing a plan view of the open well trajectory from said 

15 directional survey data, selected at least three points on the trajectory as centres of circles with the 15 

radius of each circle equal to the previously determined distance of its centre from the cased well and 
determining the location of the cased well relative to the open well by the common intersection of said 

circles* 

1 1 A method for determining distances and directions from an open well to a cased well 

20 comprising determining the resistivities and thicknesses of the earth strata penetrated by said well, by 20 
induction logs run in said open well or previously run in said cased well, determining the actual 
trajectory in the earth of the open well by directionally surveying said well, determining theoretically 
expected resistivity curves for selected resistivity tools at various depths in said open well for different 
assumed positions of the cased well relative to the open well, from said earth strata resistivities and 

25 thicknesses determined by induction logs, determining measured resistivity curves for said selected 25 
resistivity curve matches the measured resistivity curve in a least squares fit mathematical sense, and 
determining the distances of selected points in the open well to the cased well from the amount of 
translation of said theoretically expected resistivity curve required to produce said match of said 
theoretically expected and said measured resistivity curves. 

30 12. A method according to any one of claims 1 to 1 1 in which said directional survey data are 30 

obtained during drilling. 

1 3. A method for determining distance and directions of a cased well from an open well 
according to claim 1 substantially as described hereinbefore with particular reference to the 
accompanying drawings. 
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